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COMPUTER PROGRAMS FOR THE INTERPRETATION OF LOW RESOLUTION 


MASS SPECTRA: PROGRAM FOR CALCULATION OF MOLECULAR 
ISOTOPIC DISTRIBUTION AND PROGRAM FOR ASSIGNMENT 
OF MOLECULAR FORMULAS 
by Robert A, Miller* and Fred J. Kohl 
Lewis Research Center 

SUMMARY 

Two FORTRAN computer programs for the interpretation of low resolution mass 
spectra have been prepared and tested. One is for the calculation of the molecular iso- 
topic distribution of any species from stored elemental distributions. The program re- 
quires only the input of the molecular formula and has been designed to be compatable 
wife any computer system. The other program is for the determination of all possible 
combinations of atoms (and radicals) which may form an ion having a particular integer 
mass. It also uses a simplified input scheme and has been designed to be compatable 
with any system. 


INTRODUCTION 

The mass spectroscopist is faced with several laborious tasks when he is attempt- 
ing to identify an ion from its mass and isotopic distribution. One task is to select all 
ions which may give rise to a particular mass. Another is to calculate the relative 
molecular isotopic abundances of these ions for comparison with fee experimental dis- 
tribution. Two new computer programs which may be used for these tasks are de- 
scribed here. These programs are easier to handle and more encompassing than pre- 
vious methods and computer programs (refs. 1 and 2). The ease of use of our pro- 
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grams arises from having all necessary atomic data contained within the program. 
Moreover, the programs are easily adaptable to most computer systems. 


CALCULATION OF MOLECULAR ISOTOPIC DISTRIBUTIONS 

Grimley (ref. 3) has presented an algorithm which may be used to obtain molecular 
isotopic distributions from known elemental distributions. With this method, two ele- 
mental distributions are first combined. The resulting diatomic distribution is then 
combined with another elemental distribution to construct a triatomic distribution. This 
process is repeated until the distribution for the desired ionic species is obtained. This 
method, althou^ tractable for hand- calculations of small ions, becomes very tedious 
for larger species. 

For multiatomic species high-speed computer techniques are most convenient. For 
this reason we have prepared a computer program which may be used to calculate these 
distributions. In order to encourage the use of this program over hand calculations, 
thereby avoiding the associated error and tedium, a simplified input scheme is used. 

The only input required is the molecular formula; all other information is stored per- 
manently within the program. For example, the distribution for the molybdenum tri- 
oxide trimer molecule (MoOg)g, would be obtained after entering MO 3 0 9. 

The output of this program is shown in figure 1(a). It consists of the molecular 
formula, the individual elemental distributions, and the molecular distribution. This 
molecular distribution is displayed both numerically and graphically; the latter display 
is especially useful for rapid comparisons of experimental and calculated distributions. 
The run time for this calculation was 1.2 seconds on the IBM 360. A second example, 
for dibromochlorotrifluorobutane, C^H^BrgClFg, is given in figure 1(b). 

The few other programs available for this type of calculation (refs. 1 and 2) all 
have limitations. Of the three programs given in these references, one is limited to a 
few atoms; another requires the input of atomic distributions for each atom in the 
molecule; and the remaining program requires the input of a complete table of distri- 
butions with each run. Furthermore, all require the input of additional data such as the 
atomic numbers and the number of atoms. 

A listing of our program is given in appendix A. This program has run success- 
fully on the UNIVAC 1110 system at Lewis Research Center and, with a few modifica- 
tions, on the IBM 360 system. This program is called DISMIS for dis tribution of 
molecular topes. It is coded in standard FORTRAN IV and can be readily adapted to 
both time- sharing and batch-processing systems. Moreover, its output can be easily 
changed from the 129-column form to a 79-column form for teletypes. The required 
modifications are given in comment cards located within program listing. The ele- 
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mental isotopic distributions of all stable elements are contained in a subroutine called 
DISEIS. These data were taken from a recent tabulation by Roboz (ref. 4). 

The program is set up to accept molecules having up to 10 different atomic species. 
The molecular formula, which is the only input required, is entered in 10(A2,I2) for- 
mat. The common symbols for the elements are used, and those symbols and sub- 
scripts having a single character must be preceded by a blank space. Sufficient mem- 
ory has been set aside within the program to calculate the distribution for virtually any 
molecule which can be detected mass spectrometrically. 


ASSIGNMENT OF MOLECULAR FORMULAS 

The technique described here can be used to prepare a list of those species, pos- 
sibly present in a system, which have a specified mass number. For example, in fig- 
ure 2(a) all possible combinations of K, Cl, Cr, and O atoms which have a mass num- 
ber of 294 have been determined. From this list one may choose those species which 
are intuitively reasonable, and then further analyze them using the isotopic distribu- 
tion program. For the example in figure 2(a) the best candidate was the potassium di- 
chromate ion, K2Cr20^^. The run time required for this calculation was 1.2 seconds 
on the IBM 360. A second example is shown in figure 2(b), where the most reasonable 
choice for the mass to charge ratio m/e 118 was CrO 2 ( 0 H) 2 ^- 

There have been a few tabulations of molecular formulas as a function of mass 
(refs. 5 and 6), but these are limited to four and 16 elements, respectively. There are 
also a few computer programs available (e.g. , refs. 7 and 8) . However, of the three 
programs listed in these references, two are designed for high resolution, generally 
organic, spectra, and the third is written in BASIC. All are limited to a few atoms. 

A listing of our program, including pertinent COMMENT cards for the conversion 
from batch processing to time sharing computer systems, is given in appendix B. This 
version has run on the IBM 360. The mass number of the most abundant isotope for 
57 elements have been included in the program. These represent all of the stable ele- 
ments for which this isotope accounts for at least 60 percent of the total. Thus, ele- 
ments such as platinum have been excluded. In addition to the elements contained in 
this program, the radicals OH, H 2 O, and CHg have also been included. Up to 12 more 
radicals may be added to the array ATOM by replacing two of the present zeroes with a 
two-digit symbol and the mass number, respectively. More radicals may be added if 
the dimension is increased. 

The input to the program consists of the symbols of up to nine of the stored atoms 
or radicals entered in 9A2 format. This is followed by the value of m/e under con- 
sideration. It is entered in 14 format and must be right justified; for example, a two- 
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digit number must be preceded by two blank spaces. One then enters the maximum 
number of each atom to be considered. This is mandatory with the batch processing 
program; however, with the time sharing program, upper limits for these values will 
be calculated and then lower values may be added if desired. 

As seen before, the output of the program consists of a listing of all possible com- 
binations of the entered atoms which have a specified mass. Because these molecular 
masses are based on the masses of the most abundant atomic isotope, errors can occur 
if those atomic isotopes do not account for most of the intensity. This is especially 
true when several of the same atoms are being considered. For example, for the ion 
K^Clg , the atoms K and Cl account for about 93 and 75 percent of the total, re- 
spectively. For this ion the most abundant mass is not 261 atomic mass units (amu) as 
one may have expected but instead 263 amu. Thus, one must always consider the full 
distribution, as calculated using the program DISMIS, before assigning an ion to a mass 
number. 

The user of this program may also wish to limit both the number of different atoms 
to be considered and the maximum number of each. Such restrictions, when made 
carefully, can often result in substantial savings of computer time without missing any 
reasonable species. However, the user must severely restrict the number of hydrogen 
atoms to be considered. If this is not done, the number of species having the molecular 
mass of interest would be prohibitive. Thus, the program is especially useful for in- 
organic species where the number of hydrogen atoms in a molecule is generally less 
than for organics. Indeed, we have used the program to interpret the low resolution 
spectra of inorganic systems and have found it to be an invaluable aid for recognizing 
all possible species which may be in a system. 

The program is called FINDME which stands for: find m/ e. As with the pre- 
viously described program it has been designed for ease of use, compatability with 
other FORTRAN IV compilers, use with batch processing or time sharing systems, and 
use with any output device having, in this case, at least 52 columns. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 28, 1977, 

506-16. 
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APPENDIX A 


LISTING OF COMPUTER PROGRAM FOR THE CALCULATION OF 


MOLECULAR ISOTOPIC DISTRIBUTIONS 


OI»1ENSION AT(10),NAT(1Q),ATM(100) 

C FOR TTY CHANGE TO PL0T(5f)) 

DIMENSION A1(500),A2U0),FRC(50D),PL0TC1DD) 

DATA DASH , BLANK , PLUS/ 

C FOR TIME SHARING SYSTEMS REMOVE' THE FOLLOWING *C* 

C URITE(6,5nO) 

500 FORMAT ( ’OENTER MOLECULE IN 10(A2,I2) FORMAT; RIGHT JUSTIFY*) 
50 REA0(5,10O»EN0:51) ( A T ( I ) , N A T ( I ) , 1 = 1 , 1 0 ) 

100 F0RMAT(10(A2,I2)) 

NlS=in 

NA = 0 

DO 1 I : 1 , N I S 
J = NAT ( I ) 

IF(J.EQ.u) GOTO 1 
I T0P=I 
DO -2 Ilri , J 
NA:NA*1 
ATM(NA)=AT(I) 

2 CONTINUE 

1 CONTINUE 

WRITE (6, 101) (AT(I),NAT(I),I=1,IT0P) 

101 F0RMAT(///*1 ISOTOPIC DISTRIBUTION FOR *10(A2,I2)) 
WR1TF<6,102) 

102 F0RKAT(*0 ELEMENTAL DISTRIBUTIONS’) 

CALL DISEIS ( ATM ( 1 ) ,L 1 , A1 , SUM,NSI ) 

LH=L1 +NSI -1 

LE = LH 

IF(NSI.GT.7) LErLl+6 
WRITE(6,103) ATM( 1 ) , ( M, M=L 1 ,LE ) 

NIS0=NSI 

IF(NSI.GT.7) MS0 = 7 
WRITE(6,107) ( A1 ( I ) , m ,NIS0 ) 

IF(NSI.LE.7) GOTO 20 
LE=LE*1 

WRITE(6,108) (M,M=LE,LH) 

WRITE(6,1D7) (Alt!) ,I:8,NSI) 

20 IF ( ( SUM.GT. 1. 0005 ) .OR. ( SUM.LT . .9995) ) WRITE(6,109) SUM 

103 FORMAT( A7,I3,6I11 ) 

107 F0RMAT<7X,7F1 1 .6) 

108 FORMAT ( 7X , 13, 61 1 1 ) 

109 F0RMAT(9X,* NOTE SUM =*F9.6) 

NSI=NIS 

ZEROrl .E-8 
DO 3 IA=2,NA 

CALL DISEIS(ATM(IA),L2,A2,SUM,NIS) 

IF( ATM( lA ) .EQ . ATM( 1 ) ) GOTO 21 
LH = L2-*-NIS-l 
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LEtLH 

IF(MS.GT.7) LE-L2 + 6 
WRlTE(6,in3) ATM(IA) , (M,M:L2,LE ) 

MISOrMS 

IF(N1S.GT.7) NIS0:7 
WRITE(6,in7) (A2(I) ,1=1 , NISO) 

IF(N1S.LE.7) GOTO 22 
LE=LE+1 

URlTE(6,in8) (M,M=LE,LH) 

WRITE(6,1G7) (A?(I) ,I=8,NIS) 

22 IF ( ISUM.GT. l.QOni ) .OR. ( SUM.LT. .9999) ) WRITE(6,109) SUM 

21 CONTINUE 

ATM(l)rATMdfl) 

L1=L1+L2 
NI=NSI+NIS-1 
DO 4 NMin,NI 
N=NMI ♦ 1 
NMT=N«I 

IFINMT.C-I.NSI) NMT-NSI 
N S = N - M S 

IF ( NS .LE .C) NS:l 
FRAC=D. 

DO S I1:NS,NMT 
I 2:N-I 1 

FRAC^Al (II )*A2(I2)*FRAC 

5 CONTINUE 

FRC ( NMI ) :FRAC 

IF(FRAC.LT.ZERO) FRC(NMI)rO. 

IF(FRAC.GT.ZERO) NSS=NMI 
4 CONTINUE 

NSI^NSS 
DO 6 11 = 1, NSI 
A 1( II ) =FRC ( II ) 

6 CONTINUE 

3 CONTINUE 

NIMIN=n 
NIMAX=NSI 
S U M ,M r a . 

PRCMAX=D. 

DO 7 I = 1 , N S I 
P=FRC ( I ) 

IF ( F . GT . FRCMA X ) FRCMAXrF 
IFCG.^.Zi’RO) NTMAXn 
SUMM = F + SU''M 
IF ( NImIN.GT.D ) GOTO 7 
IF(F.GT.ZCRO) NIMIN=I 

7 CONTINUE 

C FOR TTY CHANGE TO NPL0T=50 

NPLOT: 100 
IPLT=NPL0T/1G 
WRI T£ ( 6 , 104 ) (1,1 = 1,10) 

104 FORMAT(’D molecular DISTRIBUTION’,/,- 
C FOR TTY CHANGE TO 1015 

£ • M/F ABSOLUTE RELATIVE ’,10110) 

DO 8 I=NIMIN, NIwaX 
F=FRC ( I ) 

MASS = L1-’I-1 
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RtLF(?C = F/rRCMAX 

IF( RFLFRC.LT. .005 ) GOTO 19 

D09J-1,NPL0T 

9 PLOT ( J) =SLANK 
M:RELFRC*NPL0T+ .5 
DC 10 Jri.N 
OL0T( JlrDASH 

10 CONTINUE 

DO 11 J=IPLT,N,TPLT 
IF(J.LE.N) PLOT(J):PLUS 

11 CONTINUE 

HRITE(6,inS) MASS, F.RELFRC, PLOT 
GOTO 8 

19 WRITE(6,in5) MASS,F,RELFRC 

8 CONTINUE 

C POR TTY CHANGE TO 50A1 

105 F0RMAT(I6,F11.7,F11.7,2X,100A1) 
WRITE(6,106) SUMM 

106 F0RMAT(*0 TOTAL PROBABILITY =*,F9.6,///) 
GOTO 50 

51 STOP 

END 


SUBROUTINE DISEIS(AT,L,A,SUM,NIST» 

DIMENSION A(500),AT0M(15,85),ATAM(15,19),ATEM(15,19) 
dimension ATIM(15,15»,ATUM(15,18),ATYM(15,14) 
data ATAM/ 

C ’ H' , 1 ., .999355, .000145, 1 1*0. , 

C ’HE* ,3., .00000137, .99999863,11*0. , 

C 'Ll ’ ,6. , .0750, .9250, 11 *0 . , 

L ’BE ’ ,9. , 1 . , 12*n. , 
t • B ’, 10. ,. 1976 , .6022 , 1 1*0. , 
t ’ C’ , 12. , .98888, .0 11 12 ,1 1*0. , 
t ’ N’ , 14 ., .99633, .00367 , 1 1*0. , 

C ’ 0’ ,16. ,.99759, .000374, .002039,10*0., 

E ’ F • , 19. , 1 . , 12*0. , 

C ’NE’, 20. ,.9092, .00 257, .0 882,10*0. , 

E ’NA’ ,23. ,1 . , 12*0. , 

E ’MG’, 24. ,.7870, .1013, .1117,10*0., 

E • AL ' ,27. , 1 . , 12*0. , 

E ’SI’, 28. ,.92 21, .0470, .03 09, 10*0., 

E • P’,31.,1., 12*0. , 

E ’ S • ,32. , .95018, .00760 ,. 04215,0. , .00014, 3*0. , 

E ’CL’ ,35. ,.7553,. 0, .2447,10*0., 

E ’AR’, 36. ,.00337,0. ,.00063,0. ,.99 600, 8*0., 

E • K’ ,39. ,.9310, .000118, .0688, 10*0./ 

DATA ATEM/ 

E ’CA’, 40. ,.9697,0. ,.0064, .00145, .020 6,0. ,.00003,0. ,.001 85, 4*0., 
E ’ SC’ ,45. , 1 . , 12*0. , 

£ ’TI’, 46. ,.0793, .0728, .7394, .0551, .05 34, 8*0., 

E ’ V’ ,50. , .0024 , .9976 ,1 1*0. , 

E ’Cfi’, 50. ,.0431,0. ,.8376, .0955, .02 38, 8*0., 
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£ ‘MN* ,55. , 1 . , 12*0. , 

t •rE*,54.,.05 82,0.,.916fc,.021<?,.0033,fl*0., 

C *C0‘ ,59. , 1 . , 12*CJ. , 

L *NI* ,58., .6788, 0.,.2623,.Q119j. 0366,0., .0108, 6^u., 
t *CU’, 63. ,.69 09,0. ,.3091, 10*0., 

t *ZN*, 64. ,.48 89,0. ,.2 78 1, .04 11, .1857,0. , .0062,6*0. , 

€ *GA* ,69. , .604 ,n. , . 396, 10*0. , 

£ *5E* ,70. , .2052 ,0. , .2743, .0776, .3654,0. ,.0776, 6*0. , 

£ 'AS* ,75., 1., 12*0., 

£ *SE ’ ,74 . , .0087,0. , .090?, .0758, ,2 352,0. , .4982,0. , .0919 ,4*0., 

£ *BR', 79.,. 50537,0.,. 49463, 10*0., 

£ *KR*, 7 8., .00 35,0. ,.0227,0. ,.11 56, .11 5 5, .5690,0. ,.1737, 4*0., 

£ ’RE *,85. ,.7215,0., .2785,10*0. , 

£ ‘SR*, 84. ,.0056,0. ,.0986, .070?, .8256, 8*0./ 
data ATIM/ 

£ • Y • , 89. , 1 . ,12*0. , 

£ ‘ZR*, 90. ,.5146, .1123, .1711,0. ,.1740,0. ,.0280, 6*0., 

£ ‘NB* ,93. ,1 . , 12*0. , 

£ ’MO’, 92. ,.15 84,0. ,.3904, .157?, .1653, . 094 6 , . 2 3 78 , 0 . , . 0 96 3 , 4*0 . , 

£ ’RU’, 96. ,.0551,0. ,.3187, .1272, .1262, .1707, .3163,0. ,.1858, 4*0., 

£ 'RH* ,103. ,1 . ,12*0. , 

£ 'PD*, 102. ,.0096,0. ,.1097 , .222 3, . 2733,0. ,.2671,0., .1181, 4*0., 

£ ’AG’, 107. ,.51817,0. ,.48183, 10*0., 

£ 'CD*, 106. ,.0122,0. ,.0088,0. ,.1239, .1275, .2407, .1226, .2886,0., 

£ .075 3,0. ,0., 

£ 'IN‘, 113. ,.0428,0. ,.9572, 10*0., 

£ *SN’, 112. ,.0096,0. ,.0066, .0035, .1430, .0761, .2403, .0858, .328 5,0., 

£ .0472,0. , .0594 , 

£ »SB’, 121. ,.5725,0. ,.42 75, 10*0., 

£ *TE’, 120. ,.00089,0., .0246, .0087, .0461, .0699, .1871,0. ,.31 79,0., 

£ .3448,0. ,0., 

£ • !• ,1?7.,1.,1?*0., 

£ •XE’,1?4.,. 00096,0. ,.00 090,0. ,.01919, .2644, .0406, .2118, .2689, 

£ 0 . , . 1044 ,0. , .0887/ 

data atuk/ 

£ •CS’,133.,1. ,12*0. , 

£ ’9A’, 130. ,.00101,0., .0 009 7,0. , . 0242, . 06 59, .0781, .1132, .71 66, 4 

£ 'LA*, 138. ,.00089, .99911, 11*0., 

£ *CE’ ,136. ,. 00103 , 0 . ,.00250,0. ,.83 48,0. ,.1107, 6*0., 

£ ’PR', 141. ,1. , 1?*G. , 

£ • ND’, 142. ,. 2711, . 1217, . 2385, . 08 30, . 172?, 0.,. 0573,0. ,.0562,4*0., 

£ ’SM’, 144. ,.0309,0. ,3. ,.1497, .1124, .1383, .0744,0. ,.2672,0., 

£ .2271,0. ,0., 

£ *EU' ,1 51., .4782,0. ,.5218, 10*0., 

£ ’GD', 152. ,.0020,0. ,.0215, .1473, .2047, .156 8, .2467,0. ,.2190, 4*0., 

£ ’TB* ,159. ,1. ,12*0. , 

£ *DY’, 156. ,.00052,0. ,.00090,0. ,.0229, .1888, .2553, .2497, .2818, 4*0., 

£ ’HC* , 165 . , 1 . , 1?*0. , 

£ 'ER', 162. ,.001 36,0. ,.0156,0. ,.3341, .2294 , .2 70 7,0. ,.1488, 4*0 ., 

£ *TK* ,169. ,1. ,12*0. , 

£ •YB’tl66. ,. 00135, 0.,. 0303, .1431,. 2182, .1613, .3164,0. ,. 1273, 4*0., 

£ ’LO* , 175 . , .9741 , .0259, 1 1*0., 

£ -HF', 174. ,.0018,0. ,.05 20, .1850, .2714, .1375, .3 524, 6*0., 

£ • T A 1 80 ., .00012 , .99988, 1 1*0./ 
d/ita ATYM/ 

£ ’ W’, 180. ,.00 14,0. ,.2641, .1440, .3064,0. ,.2641, 6*0., 

£ 'RE*, 185. ,.370 7,0. ,.6293, 10*0. , 


8 



C *0S* flB4* f •OGG?fO« f*0159f •0164t* 133t • 161f *2b4f0« y y^lGy^^Oyf 
L yl91 .373,0. , .627,10*0. t 

t ’PT*, 190. ,.000 13,0., .0078,0. ,. 329, .338, .253,0. ,.0721, 4*0., 

C ’AU’ ,197. ,1 . ,12*0. , 

£ 'HG’, 196. ,.00146,0. ,.1002, .1684, .2313, .1322, .2980,0. ,.0685, 4*0., 
t *TL’ ,203. , .2950,0. , .7050,10*0. , 

£ ’PB* ,204 . ,.0148,0. ,.236, .226, .523, 8*0., 

£ *31’ ,209., 1 . ,12*0., 

£ ’TH* ,232.,!. ,12*0., 

£ ’ U’ ,234 ., .000057 , .0072, Q. ,0. , .9927, 8*0. , 

£ 30*0./ 

E QUlVALENCE(ATOM(l),ATAM(in,< ATOM <286), ATEM(1)),<AT0M(671), 

£ ATIM(1)),(AT0M(796),ATUK(1M, < ATOM ( 1066) , ATYM ( I ) ) 

SUM=0. 

ATM=AT 
DO 1 1=1 , 85 
J = I 

ir ( ATM.EQ . atom ( 1 , I ) ) GOTO 2 

1 CONTINUE 

2 L=IFIX(AT0M(2,J)) 

0031=1,13 

A(I)=AT0M(I+2,J) 

IF < A ( I ) .LE . 0. ) GOTO 3 
SUM=SUM+A(I) 

NIST=I 

3 - CONTINUE 

PETUPN 

END 
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APPENDIX B 


LISTING OF COMPUTER PROGRAM FOR THE ASSIGNMENT 
OF MOLECULAR FORMULAS 


INTEGER AT0M(1*40),AMASS(9),ATM(9),L(9),J(9) 

DATA ATOM/* H * , 1 , * L I * , 7 , * 9 E * , 9 , • B*,ll,* C*,12,* N * , 1 4 , * ME ’ , 1 5 , 
t * 0* ,16,’OH* ,17,*UA*,ie,* F * ,19, *NE* ,20, *NA • ,23, ’MG’ ,24, ’ AL * ,27, 
L *SI*,2e,* P’,31,* S* , 32, *CL* ,35, * AP* ,40, * ,K* ,39, *CA * ,40, * SC * ,45 , 
£ *TI*,49,* V*,51,*CR*,52,*MN*,55,*FE’,56,*C0*,59,*NI*,58,*CU*,63, 
L *GA* ,69, * AS* t75,* RB* ,72, *SR* ,88, * Y * , 89 , * NB * , 9 3 , * RH * , 1 03 , 

£ *IN*,96,* I*, 10Q-, *CS* , 133,*BA*,138,*LA*,139,*CE*,140,*PR*,141, 

£ *TB*,159,*H0*,165,*TM*,169,*LU*,175,*TA*,191,*RE*,187,*IR*,193, 

£ * AU* , 197 , *TL * ,205, *31 * ,209, *TH*, 232, * U*,238, 

£ 24«0, 

£ * *, 0 / 

DATA YES/*Y*/ 

C ^OR TIME SHARING SYSTEMS REMOVE *C* FROM I/O AND IF STATEMENTS 

C WRITE(6,200) 

200 FORMAT ( ’CENTER ATOMS AS 9A2*) 

501 READ(5 , 10G,E.N0:5QG ) ATM 

100 F0RMAT(9A2) 

C WRITE<6,2ni) 

201 FORMAT( ’CENTER M/E AS 14*) 

READ(S,101) MOLMAS 

101 F0RMAT(I4) 

00 1 1=1,9 

DC 2 K=l,99,2 

IF < ATM( T ) .NE . AT0M(K n GOTO 2 
A MA SS ( I ) =A TOM ( M -*1 ) 

GOTO 1 

2 CONTINUE 

1 CONTINUE 

102 FORMAT ( ’DMASS *,9A3) 

DO 5 1=1,9 

NB=I 

IF ( ATM( I ) .FQ, AT0M<99) ) GOTO 6 

5 CONTINUE 
NB = N3* 1 

6 CONTINUE 

17 NC=N3-1 

NCMR=1 

DO 16 IT=1,NC 
L (II)=MOLMAS/AMASS(II ) 

NCMB= (L ( 1 1 ) ♦! )«NCM3 
16 CONTINUE 

WRITE(6,202) NC MB , ( L ( I ) , I = 1 , NC ) 

202 FORMAK* THERE ARE’.IID,* COMBINATIONS*,/ 

£ * MAX NO. OF EACH ATOM IS*, 913) 

C JRITE (6 , 109 ) 

109 FORMAK* DO YOU WISH TO RESET MAXJMUMS *) 

C REA0(5,2C3) HUH 
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203 FORMAT(ftl) 

C IF(HUH.Nl.YES) goto 24 

C WRITE(6,204) 

204 FORMAT( ’CENTER NO. ATOMS AS 912’) 
READ(5,20S)(L(I),I=1,NC) 

205 F0RMAT(9I?» 

24 IMASS=n 

L1=L( 1 )+l 
L2 = L(2)-’l 
L3 = L< 3) ♦! 

L4 = L ( 4 )4l 
L5 = L( 5)-’l 
L6=L(6)+1 
L7=L(7)+1 
L8rL(8)*l 
L9=L(9)+1 

WRITE«6,102)(ATM<I),I=1,NC) 

DO 30 11=1, LI 

J(1 )=I1-1 

41 = J( 1 )*AMASS(1 ) 

DO 30 12=1, L2 
J(2)=I2-1 

M2=J( ?)OAMASS (2)+Ml 
00 30 13=1, L3 
J ( 3 ) =1 3-1 

M3 = J( 3)*AMASS(3)-’M2 
DO 30 14=1, L4 
J(4)=I4-1 

M4 = J< 4)*AMASS (4) *M3 
00 31 15 = 1 , L5 
J<5»=15-1 

M5=J( SlOAMASS (5)^H4 
IF ( M5 .GT. MOLM AS ) GOTO 30 
DO 32 16=1, L6 
J<6)=16-1 

M6 = J( 6 ) *AMASS < 6 ) •’MS 
IF(M6.GT.M0LMAS) GOTO 31 
DO 33 17=1, L7 
J(7)=I7-1 

M7 = J( 7 >*AMASS ( 7 »-’M6 
IF(M7.GT,M0LMAS) GOTO 32 
DO 34 18=1, L8 
J(8)=I8-1 

M8 = J(8)*AMASS(8)’-M7 
IF ( M8 .GT. MOLMAS ) GOTO 33 
DO 35 19=1, L9 
J (9 )=T9-1 

MSS=J(9)>«‘AMASS(9)+M8 
IF( MOLMAS-MSS) 34,21,35 
21 WRITE(6,108) MS S , ( J ( I) , I = 1 , NC > 

108 FORMAT! 15 ,913 ) 

35 CONTINUE 

34 CONTINUE 

33 CONTINUE 

32 CONTINUE 

31 CONTINUE 

30 CONTINUE 



C PRINT OUT OF MAX. NO. OF EACH ATOM (LEAVE OUT OF T/S VERSION) 

WRITE(6,104) (L(l) ,1=1,NC) 

104 F0RMAT(/5X,9I3,//) 

SOTO 501 
50Q STOP 
END 
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ENTER MOLECULE IN lnCA2,l2) FORMAT; RICIIT JIICTIFY 
mo j o 3 


♦ ♦ 


h- rn O o o o 

•j>cr> . o . .-4 . 

C3 o o o 


<t rj i-i Cl o o cr» 

f-" CD *0 • rH • 

z: o o C3 


LJ lA un J «H O O C-4 to to O :-•% O O ^ rH C3 O C5 O iH O lA r.*l CJ rj 
> O ;A ^ O C*. f-4 to UO O !>• CO O lA rA lA i-t r-t r«. t*l O 50 O U*\ hA O r-4 

— ocNr^r-^tACMrvtoc^o^c: ooots*c 4 >»' 3 r>.rAr^cocirjrAOOO 

H «-» r-» o »A O if O O O 3 <4 r-l t-O -A O r** O O U> O r-i O O 

z: < r>. O ^ O ^ O CO C4 O -f LA c*l CO O fA r-t O .“J fA tO --fl '.0 O O D O O 

c — J m 3 o r-« .A rs C'4 .*■1 c» w o »A -0 in A CO rN. o ^4 o 3 3 o o o o 

— LJOO 3 r-l«Hf-»^»A:- 0 'wrN.t 0 Or^C 0 t 0 O-^IArHf' 4 OOOOOf 3 OO 

zv 30 0 00300 otoocor-toooocos CO oooroooaoo 

CC UJ C3 fH Cl C>l >1 O c^4 CM JA f-l c-o J- O rA O O C4 •/•I to 1-1 CO if if O CM rv. I-I O iH 

^-C3rnOi— iT'.'OCTf-i 3 i-i if ?s. tA if i-l o o -A C) CO O o : T o r A >• o o 

•.o _Z to «-H o • A lA CO t-4 lA r«. o ^ r-i CO JA o O to c> o o r>.» o o tH o o 

— -JC.0O'*«.Oi-«ifO-3'OO 3r-«:AtAOOrArHr>.CACN.CICOOOOO oo 
C C fA O O .H o O »A O -A O lA .f if -.0 CO O CO r«» r-l r>. tH .-Nl LO O O O O O O 

oo o o*H«-«tMifrAors.to^oor^orN.oif-fiHf4 0 o oooooo 

ZCC^OOOOOOOOOOOOr-IOOOOOOOOOOOOOOOO 

< 

_l OOOOOOO 0-0 000 0 0 000000000000000 

3 

o 

LJ LJ -O r-4 Al fA if '-A O r«. CO O O rH *4 }A if tA lO r» CO Cl O iH AJ r-A if lA O r«.. -0 
_l ■>.. CM Al Al CM CM *N| CM .M T-l CM :A rA fA rA lA fA fA lAt-AtAif-f-Cfif-fifif^J- 
OlIififif-f^if-f-f-fif-fifififJ-ifififif«f-^’-0--fifififififif 
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(a) Molybdenum trloxide trimer. 



EIITER MOLECULE iri 10(A2,I2) FORMAT; RIOIIT JUSTIFY 


O 


^-^“»o^^ooooo 

ZD . 


UJ-^coofOrHinrvsor^cD 

>«NCNO»>.r*.OCMCTLDO 

— u>oo^cr»coKvfr»oo 
I— »-ioOLni-'\r«.mf^i— to 
z:<cocr>omcoi-toooo 

0 ~-K\i- 40 ^CTJ^rOOOO 
— UJ^OOOOOrHOOO 

HC 

13 OOrHOOOOOOO 


CSMCTf-IC»C>lAtAr».00 

• rH • • • 


•f— rs.C»f- 4 J-»-«lAUi 3 JS '0 
> 3 >tf\cr»'-r\omr^J’ o^o 
- _l » A r>. rH OO - 3 - O O o 

iO-^caoOtOTAr^rjoo 
CO»OCMrHOrHinOOO 
:CD«-HOJ-OCNJOOOOO 

c < 

J OOOOOOO 300 


I lU rH K\ ^ LP ^ O O rH 
J-VsOO'OOOOOOl— tiH 
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run f i ndne 


ENTER ATOMS AS 3A2 
kc 1 c r o 

ENTER H/E AS 14 
294 

THERE ARE 8208 CONCIMATIONS 

MAX NO. OF EACH ATOM IS 7 8 5 18 

DO YOU WISH TO RESET MAXIMUMS? 

no 

MASS K CL CR 0 

294 U 2 0 14 

294 0 2 4 1 

294 0 0 1 2 

294 1 1 3 4 

294 1 5 l) 5 

294 2 U 2 7 

294 4 2 1 1 

294 5 1 0 4 


(a) Potassium, chlorine, chromium, and 
oxygen for mass of 294 atomic mass 
units. 


run f Inline 

ENTER ATiirS AS 0A2 

cr onac 1 oha 1 s I fen I 

ENTER N/E AS 14 
118 

THERE ARE 907200 CONC I NAT I ONS 
MAX NO. OF EACH ATOI' I S 2 7 5 3 C 4 4 2 2 

00 YOU WISH TO RESET MAXIMUMS? 

yes 

ENTER NO. ATOMS AS 9 1 2 
2 7 4342211 

M.ASS CR 0 NA CL OH Al. SI FE III 
llo 000020110 
118 000101010 
118 000101200 
118UU1041000 
118 002011100 
118 0 1 0 0 1 1 0 0 1 

118 011030100 

118 011111000 
118 012000010 
118 012000200 
118 020000101 
118 020130000 

118 0 2 1100100 
118 0 2 3010000 

118 0 3 0 2 0 0 0 0 0 

118 0 4 0 0 0 2 0 0 0 

118 1 1 1 0 0 1 0 0 0 

118 120020000 


(b) Chromium, oxygen, sodium, chlorine, hydroxyl, aluminum, 
silicon, iron, and nickel for mass of 118 atomic mass units. 


Figure 2. - Combinations of elements or radicals. Output taken 
from time-sharing version of progran^ FINDME . 
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